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■ABSTRAC T 


The isobar ic vapour-liquid equilibrium data have 
been measured for the system Methylcyclohexanol + Methylcyclohexane. 
03ie equilibrium still, as designed by Yerazunis, Plowright and 
Smola has been used for this purpose after some modifications. 

The data have been tested for thermodynamic consistency using 
the integral test. Also, a method for the prediction of iso- 
thermal vapour-liquid equilibrium data from isobaric vapour 
liquid equilibrium data has been developed. It is based on the 
introduction of temperature dependence of the constants in 
Redlich-Kister equation for excess Gibb's energy. It has been 
found satisfactory for simple as well as complex mixtures. 



CHAPTER 1 


IFTEODUCTIOH 

Experimental vapour liquid equilibrium studies still 
continue both for the purpose of collecting data useful for 
design of industrial distillation units and for a theoretical 
understanding of the behaviour of liquid mixtures. 

The present investigation deals with the vapour liquid 
equilibria of hydrocarbon + alcohol mixtures. According to the 
work done so far it is generally found that the azeotropes are 
formed in such mixtures when there is a large difference in 
molecular size and structure of the two components. Examples 

1 r 

are Methanol + Toluene- (azeotrope at ,88 mole ^ MeOH) , Ethanol- 

1 7 1 R 

Benzene (azeotrope at 50 mole ^) , Propanol-Benzene (azeotrope at 

35 mole ^ PrOH), Propanol-Cyclohexane^ (azeotrope at 60 mole fo PrOH) 
Propanol-Hexane (azeotrope at 5 mole io PrOH) etc. However, if 
one wishes to isolate the effect of H - OH interactions on the 
vapour liquid equilibrium data one should study mixtures of compo- 
nents having same level of molecular size and structure. Studies 

carried out in this direction in the literature are for systems 

21 21 2 7 

Toluene + Phenol, Me thylcyclohexane + Phenol, Benzene + Cresol.'" 

23 ’ 

Recently Jffukherjee studied the following systems Methylcyclohexano 
Heptane and Eeptanol-Toluene3 These do not shov.’ any azeotropic 
behaviour. This work is essentially a continuation in this 
direction. The system chosen is Methylcyc lohexanol+Methy Icyclo- 
hexane. Prom the behaviour pf the above systems it is expected 



2 


that the system Methylcyclohexanol + Methylcyclohexane will not 
show azeotropic behaviour. 

Isobaric vapour liquid equilibrium data have been measured 
for the above mentioned system. The equilibrium still used for 

5 

such determination is the one designed recently by Yerazunis, 
Plo¥/right and Smola. This works on circulation of both* vapour 
and liquid phases. The highlights and special features of this 
still are given in the second chapter. 

Isobaric vapour liquid equilibria does not directly give 
a clear picture of non-ideality and behaviour of excess G-ibbs 
energy of mixing. The information is rather diffused. However, 
a method has been developed to predict isothermal data from 
isobaric data and is based on the introduction of temperature 
dependence of the constants in the Redlich-Kister equation for 
excess Gibbs energy. The predictions agree well with experimental 
data for a number of systems, and in this process it is revealed 
that isobaric data if handled carefully can give lot of informat io 
regarding the behaviour of non-ideal solutions. 





CHAPTER 2 


EXPERIMTAL SET-UR MR TROCEEUBE; 

* 

The need for accurate vapour liquid equilibrium data has 

given rise to a large number of equilibrium stills - the apparatus 

in which vapour pressure, temperature and composition are usually 

measured. Before describing the equilibrium still used in the 

present study a brief history of the development of the equilibrium 

1 

stills will not be found out of place. It was Othmer who success- 
fully operated a still of the vapour recirculatory kind. Other 
workers made several other modifications. However, heat losses 
leading to partial condensation, inadequate mixing of the recycled 

vapour condensate with the liquid in the re boiler still remained 

2 

the major drawbacks. Jones et.al. devised another still of the 
vapour recirculation type in which heat losses in the vapour tube 
were compensated by external heating, mixing in the boiler vi&s 
promoted by bubbling the returning vapour through a nozzle in the 
re boiler. High internal pressure drop, inexact measurement of 
equilibrium temperature, and the need for exact compensation of 
heat losses (undercompensation led to stoppage of circulation and 
over compensation led to overloading of the vaporizing tube with 
the liquid) checked any further attempts at modification of 1h.e 
still in this direction. Gillespie'^ designed a still which over- 
came many of the defects discussed above. It was based on the 
vapour-re ctrculato.ry principle but differed from other types in 
that the boiling liquid is also circulated rapidly by employing 
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the Cottrell device‘s, Landwehr ,et .al . modified Gillespie 
apparatus - hy separating the point of liquid sampling from the 
re boiler, lagging the vapour In the separation chamber by recycle 
liquid stream.The present still is a highly improved version of 
the Gillespie still, overcomes certain basic flaws in the landwehr 
design also. This has been developed at Rensslaer Polytechnic 
Institute by Yerazunis, et.al. 

The Equilibrium Still (Pig.1) ; 

This remains basically a vapour and liquid recirculation 
device. Its most distinguishing feature is to be found in' 
the details of the separation chamber (Eig.2). The liquid and 
vapour are not permitted to separate within the boiler, but are 
maintained in intimate contact as they rise up the Cottrell pump 
in the form of slugs to the separation chamber. Vapour and 
liquid mixture at the top of the Cottrell pump is directed down- 
ward by a reverse bend through a short column packed with glass 
beads. The mixture gets separated only after it has passed througi. 
the packed section. The liquid streams past a thermocouple in 
the thermowell and overflows to the receiver-mixer bypassing the 
liquid sampling cell. The vapor flows upward around the packed 
column and thence downward around the entire separation chamber 
providing a thermal barrier to the separation region. The vapour 
stream then passes into a condenser, the condensate from which 
flows to the receiver-mixer bypassing the vapour condensate sampling 
cell. The stream of hot liquid disengaged from the vapor within 
the separation chamber is combined with the recycling condensate 
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before It returns to the boiler in the receiver-mixer. The 
mixing is enhanced by increasing the interfacial area by putting 
glass beads in the mixer. The use of a magnetic stirrer has 
been avoided in this way. The main features of this still ares 

( i) equilibrium temperature is measured, after ensuring 

equilibrium between the boiling liquid mixture and its 
vapour at the point where the two disengage in the 
separation chamber and not in the re boiler* This 
avoids measurement of temperature of superheated 
liquid and ensures measurement of the correct 
equilibrium temperature. 

(ii) packed column in the separation chamber provides 

additional (in addition to the Cottrell pump) contact 
for the two phases to approach equilibrium. Because 
the vapour and liquid in the Cottrell pump exist in the 
form of slugs the equilibrium is attained there only 
on a local basis. The packed column, by providing 
additional interfacial area, serves to smear out the 
local fluctuations and raises the equilibrium to an 
integral basis 

( iii) partial condensation and reflux of the vapour in the 
separation chamber is prevented by using the vapour 
phase as a thermal barrier.Iandwehr has used the 
liquid recycle stream for lagging. But if we look 
at it more closely there is a positive gain in using 
vapour for lagging. In vapour lagging latent heat of 
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condensation is involved and in liquid lagging 
the sensible heat is involved. Since latent heat 
is much more than the sensible heat, vapour lagging 
is more advantageous, 

( iv) sampling of liquid and vapour condensate is possible 
without stopping reciroulation and without disturbing 
the equilibrium. The vapor sampling cell is not 
of the conventional overflow design and therefore 
there is no enrichment of the more volatile component 
which may occur betvjeen termination of recirculation 
and sample removal. 

(v) Entrainment of vapor with liquid recycle in the 

separation chamber has been avoided by the creation 
of a liquid seal around the thermowell, 

(vi) the use of the receiver-mixer eliminates the existence.. 
of any concentration gradients in the reboiler. It 
is clear that at equilibrium the liquid-vapor mixture 
leaving the boiler is identical in quantity and 

:..5 composition with the combined liquid streams returning 

to it. 

(vii) all portions of the equilibrium still - vapour-side 
condenser, liquid-side condenser, reservoir, sample 
cells, receiver-mixer are separately connected to the 
manifold. This ensures a uniform pressure throughout 
and eliminates the chances of. internal pressure drop. 



9 


(viii) for prevention of heat Loss- to the surroundings 
the separation chamher has been surrounded by a 
double-v^alled, silvered and eyacuated jacket. 

This still differs in few aspects from the one designed 
by Yerazunis et.al. Firstly the vapour and liquid sampling loops 
have been altered. Whereas previously each loop had three three- 
way stop cocks. We have used only two ^-way stop cocks on each 
side, and the sample cells are connected to the manifold so that 
they are under system pressure. ?lien the equilibrium still is in 
operation so that the withdrawl of sample into the sample cells is 
done at the equilibrium pressure. Second alteration is in that 
the connection of the equilibrium still to the manifold is done 
through, a third additional water condenser. This prevents any 
escape of vapour into the manifold. Lastly, the drainage port 
in the reboiler has been eliminated by completely sealing off the 
reboiler. For purpose of drainage the seal is cut and liquid 
removed by means of a syringe. This change had to be made because 
the compounds under study (viz. methylcyclohexanol and methylcyclohexan' 
were eating away the grease. 

General Arrangement of the Apparatus ; (Fig. 3) 

The manifold consists of a 6 ft., I" diameter copper tube 
with six connections. Tvro connections on the left with an opening 
of 1/4" diameter are connected through a pair of toggle valves' - 
one to the nitrogen cylinder and the other to the vacuum pump. 

The third connection is via a brazed tee joint to the surge tank. 







11 


The surge tank consists of two 20 liter round bottom flasks. 

The other three connections are one each to the equilibrium 
still, manometer and ebulliometer. Equilibrium still and ebullio- 
meter connections employ ground glass B-1 2 joints so that at any 
time they can be isolated from the manifold. Vacuum tubing is 
used to make the copper-glass connections. Manometer is directly 
connected to the copper line through a vacuum tubing. 

The two 20 liter round bottomed flasks which function as 
the surge tanks are kept in a wooden box packed with insulating 
material - which lends a thermal stability to the surge tank. This 
large volume acts as a damper to pressure fluctuations that may 
arise in the system due to changes in room temperature, irregular 
heat losses, voltage fluctuations, sampling etc. 

Manometer has been used for rough measurements of pressure. 
Accurate measiirement of pressure is done by means of the single 
stage ebulliometer (Figure 4) . The ebulliometer permits exceptionali ,7 
precise measurements of boiling temperature. If we know the exact 
dependence of boiling point on pressure for a pure substance, by 
measuring the boiling points of this substance we can measure 
pressure over any range. Water, benzene, toluene are some of the 
substances for which P-T dependence is known with high precision. 

The ebulliometer (Ftg.4) essentially consists of a boiling ■ 
tube, Cottrell pump, thermowell, condenser, and a drop counter. 

The super heaiedmixture of liquid and vapour rises up the cottrell 
pump in th^^ of slugs and spurts over the thermowell which 

is wrapped with a glass spiral. The glass spiral helps in 
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breaking the liquid -vapour mixture; The liquid at this point 
looses some of its superheat and the equilibrium temperature is 
measured from the thermowell. The vapour rises up the condenser, 
condenses and recycles into the boiling tube by trickling in the 
form of drops from the drop counter. The drop rate is a measure of 
the intensity of heating. Operation is done over a range of drop 
rate, in this range the boiling point essentially remains constant. 
The lower limit signifies inadequate heating and the higher limit 
signifies excessive heating. A float trap has been provided in 
the recycle path so that the boiling liquid should not rise up on 
the other side also. 

Toggle valves are snap action valves, their opening time 
being few seconds. These valves have been used in conjunction, 
with small volume metal accumulators, they facilitate change in 
the system pressure by as small a magnitude as + ,5 mm Hg.(I’ig.5) 

Heat supply to the re boiler is done by means of two 
heaters - one external and tbe other finger heater. The finger 
heater consists of a ceramic rod over which the Nichrorae wire 
of gage 22 is wound. It has a total resistance of I 9 ohms. The 
exiternal heater is wound helically by the same nichrome wire, Th's 
provides a resistance of 25-5^ . Heating of the boiling tube of the 
ebulliometer is also done by nichrome 22 gage wire (20.^) wrapped 
helically over it. Ill the three beaters are connected to the 
main supply through variable transformers. 

^ been used for condensation. The vapour, 






Fig, 5a -Toggle Valve Stem 
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condenser, liquid side condenser and the additional condenser 
are connected in series. The sampling cell coolers form another 
set connected in series. 

Experimental Procedure ; 

Before the collection of data for the binary system 
involved (methylcyclohexanol - metliylcyclohexane) the whole set-up 
was standardized so that one could, be certain about the attainment 
of equilibrium and affect' change in the equilibrium point (say 
from x^ y^ T' to x" , y" , T") by easy control of operating conditions 
(pressure, voltage or heat input, vapor and liquid recirculation rates 
etc.). Details of standardization have been included in the 
Appendix. 

Prior to the charging of liqiid into the roboiler the whole 
system (surge tank, manifold, equilibrium still) was purged of 
air by means of a vacuum pump. When the manometer reading was 
essentially zero, suction was stopped. Inert gas nitrogen was then 
fed in through the toggle valves so that system pressure was 760 mm, 
Hg. Ebulliometer filled 77ith 60 cc of inalar Benzene was plugged to 
the manifold. The calibrated reservoir in the equilibrium still was 
filled with 200 cc of purified methylcyclohexanol. The loading port 
was closed and the two two-way stop cocks connecting the reservoir 
to the system were opened. This brought the reservoir to system 
pressure. The liquid trickles down to the reboiler along the path 
of liquid recycle stream, creating a temporary hold up in the receive 
mixer. ^-Tien all of the liquid has reached the receiver-mixer the 
reservoir is Isolated from the system. To shorten the time of total 
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transfer of the liquid to the reboiler, the heaters (external and 
finger heater both) are put on^the temperature gradient so created 
helps in quick and rapid flow to the reboiler. Simultaneous to 
the putting on of re boiler heaters, the ebulliometrtc heater is 
also put on, and the flow of water through all the condensers 
(liquid-side, vapour-side, sampling cells, ebulliometric condenser) 
started. 

Attempt is now made to achieve the equilibrium steady state. 
As revealed by standardization procedure the equilibrium is reached 
under steady condition of pressure, when (i) there is ra|d move- 
ment of slugs of liquid and vapour over the entire length of cottrell 
pump, (ii) liquid seal is formed around the thermowell- so that 
there is no entrainment of vapours with the liquid recycle stream, 
(iii) there is no hold-up in the receiver-mixer, ( iv) vapour recir- 
culation rate is high. Visual observation is made to confirm these 
facts. Voltage is adjusted to control these factors, 

for pressure adjustment recourse of IT 2 and vacuum pump is 
taken. Initially one finds that the pressure in the system increases 
This is due to expansion of 112 at the higher temperature. Manometer 
registers an increase of as much as 1 0 to I 5 mm. If the system 
pressure exceeds 760 mm, which it does normally, pressure is reduced 
b^!^ vacuum through the toggle valves. They are opened and closed 
alternately. The pair of toggle valves working in conjunction with 
the small accumulator are suitable to affect a change of +.5 mm.Hg.. 
The large volume provided by the surge tank also aids in this 
precise control. Pressure is read from the ebulliometer. At 760 mm. 
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the ebulliometer should exactly reprod.uce the normal boiling 
point of benzene, that is 80.10°C. 

When steady state regarding the system pressure has been 
reached one looks for the equilibrium steady state. The temperature 
read from the thermocouple in the separation chamber thermowell, 
should stop fluctuating, llien the potentiometer indicates a constant 
value of enif for about I 5 min. - the steady state is believed to 
have been reached. System pressure is again checked from the 
ebulliometer reading. Temperatures - equilibrium and ebulliometric 
are measured by a copper-constantan thermocouple which has heen 
previously calibrated against a platinum resistance thermometer. 

Vapour and liquid samples are now withdrawn simultaneously 
into the sample cells.. Because three-way stop-corks are used 
recirculation is not interrupted. The samples are allowed to cool 
for about five minutes. Thee ells are disconnected from the system 
by closing their connection to the manifold, and vented to the 
atmosphere. The vent prevents any contamination of the sample with 
extraneous material from outside during the time of removal. The 
samples are collected in sample tubes - the amount being approximatei 
1 cc. The samples are analysed on a ref ractometer which has been 
previously calibrated. The refractive indices measurement provides 
an accuracy of approximately to the vapour and liquid compo- 

sitions. 

The: same steady state is double checked by taking a second 
set of readings - x, y and T. Time interval kept between the two 
readings varied between 45 min, to 1 hour. 



19 


To effect a shift in composition, after the above deter- 
mination was made, a few milliliters (15 cc to 2Q cc) of the 
liauid v^ere drained through the sampling line, and roughly the same 
amount of the second component was added through the loading port. 
This procedure of changing the composition was continued till it 
v/as no longer economical i.e. the addition of the second component 
brought little or no change in the equilibrium temperature. The re- 
boiler, at this point, was emptied the system cleaned and the other 
pure component was filled in. Thus the data was collected in two 
parts - one starting from pure methylcyclohexanol and two starting 
from pure methylcyclohexane. 

About 18 data points were obtained for the entire range. 

In addition the boiling points of the pure components were also 
determined and the pure liquid and pure vapour samples for both 
the components were also analyzed on the refractometer to confirm 
the purity of the samples. 
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SOURCE MD PURITY 

Bfethylcyclohexariol was obtained as laboratorj^ grade from 
Riedel and Meth^rlcyclohexane was supplied as Analar grade from 
Eluka. Blethylcyclohexanol was purified by a process of successive 
distillations and the purit^r after each distillation was checked 
on a differential ebulliometer - this ebulliometric determination 
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(Eigure 6) of purity has been extensively used by Swe itslawsii. 

The differential ebulliometer measures the boiling point and the 
condensation temperature of the liquid in question vihen equilibrium 
conditions have been reached. The difference ( A T) between the 
two temperatures indicates the extent of purity - the smaller the 
difference purer is the compound or it has lesser impurities. 

The distillation of mothyl^yclohexanol was performed in a 
column (5 ft. column length and Internal Piometor = 40 mm) packed 
with fflass helices (average length = 1 2 mm and = 5. mm) , 

under vacuum and in an inert atmosphere of Hitrogen. The disti- 
llation head was designed for reflux - under normal conditions 
it gives a reflux ratio of 1 . However to improve the reflux 

ratio the distillate v/as collected intermittently - for 10 minutes 
it was on reflux and then for the next 10 minutes it vvas kept on 
total reflux. Details of the distillation process are being 
given below i 

Pressure of the system . = 1 00 mm - 120 mm Hg. absolute 

Temperature of the =95 to 100°C 

vapours • ■ 
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Rate of collection 

= 

20 cc/hr. 

Overall reflux ratio L/D 


2i1 

Total amount charged 

= 

1500 ml 

Total amount recovered 


800 ml 

io Recovery each time 

= 

80 to 

Rumber of distillations 
per formed 

=: 

4 

Total time spent for 
d ist illation 

= 

200 hrs. 

A - 

= 

'^b “ ^c 

T of the orginal 

rr: 

171 .19 - 167.64 

sample 

= 

3.55°C 

A T after 1st distilla- 


170.46 - 167.65 

tion 

= 

2 .81 °C 

A T after 2nd distilla- 


169.83 - 167.64 

tion 

=■ 

2 . 19 °C 

A T after 5rd distilla- 

= 

168.84 - 167.67 

tion 

= 

1 .1 7°C 

A T after 4th distilla- 

:r: 

168.42 - 167.68 

t io n 

= 

,74°C 


An essential difficulty that one faces in purifying 
methylcyclohexanol is that it forms 12 isomers (arising out of 
the combinations of orthOj meta and para positions with cis and 
trans conformations and dextro and levo rotatory behaviours). 
These isomers are spread over a large boiling range. Riedel 
boiling range is given as 165 -• 173°C. It is mot easy to 



23 


isolate these isomers just by distillation alone. Perfect 
control of pressure is also required Methylcyclohezane was also 
given an ebulliometric test. J\ T turned out to be negligible, 
less than ,1°C which is the accuracy of measurement. 

Boiling Point Boiling point Reported 

in Differential in Equilibrium boiling 

Ebul Home ter still point 


168.42 

C7H-1 ^ 


167.76 


1'65.173°C 

(Riedel) 

100.90® 


100.90 


100.71 
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EXPERmiiRTAL RESULT'S 

The isobaric experimental data obtained for the system 
methylcyclohexano 1-methylcycIohexane is reported in Table I. 

The smoothened data is reported in Table II. The accuracy of 
measurement of temperature is + .1°C and the accuracy of measure- 
ment of liquid and i7"^pour composition is ,0037. See Figure 7. 

The smoothened data has been obtained by polynomial 
regression analysis of the experimental data. The degree of 
polynomial is chooser which satisfies the experimental data 
within 95^ confidence limit and this optimum degree of the 
polynomial is obtained after applying F-testjon the experimental 
data. 

’'hen y is regressed on x the standard deviation is .006 
when T is regressed on x the standard deviation is ,22°C. The 
accuracy of measurement of composition on a refractometer as 
obtained from calibration data is .003?; and the accuracy of 
temperature measurement by means of copper constantan thermocouple, 
as indicated by the analysis of calibration data is .'1°C. ¥e 
observe that for both cases composition as- vi/ell as temperature, 
error in the experiment is greater than the accuracy. This 
indicates that there are certain inherent drawbacks in the 
experimental set-up. The experimental errors incurred will 
then affect the analysis of the experimental results viz the 
calculation of activity coefficients of both thocomponents- 
methylGyclohexanol and metlylcyclohexane . 




Fig. 7 - Vapor liquid equilibrium data tor MethylcyclohexanoKl) 

Methylcyclohexane(2) at 760mm Hg. 


26 


TABLE I 


BXPBBD, 1 ]I]ffrAL ESSUIiTS . 

at 760 mm of Eg 


X 

Y 


Pure component 

- 


.0505 

.0022 

101 .86 

.1234 

.0097 

102-38 

.1882 

.0201 

105.02 

.2556 

.0370 

104.09 

.2556 

.0343 

104.48 

.3177 

.0451 

105.82 

.5417 

.0856 

112.88 

. 5471 ' 

.0802 

112.79 

.6417 

.1072 

116.50 

.6443 

.1072 

117.02 

.7469 

.1828 

124.31 

.7496 

.1720 

124-29 

,8008 

.2287 

129.03 

.8008 

.2260 

128.81 

.9520 

.6605 

152.10 

.9601 

.7040 

154.57 

.9844 

.8629 

161 .63 

.9844 

.8629 

161 .63 

Pure component 


167.76 
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TABLE II 

SMOOTHBl'IED PJ 3 SIILTS 

at 760 mm of Hg. 


X 

Y 

T°C 

.0505 

.0025 

101 .84 

.1234 

.0096 

102 .42 

.1 622 

.0199 

103.00 

.■2556 

•0544 

104.25 

.5157 

.0492 

105.97 

.5417 

.0783 

112.68 

.5471 

.0797 

112.86 

.6416 

.1120 

116.81 

.6443 

.1132 

116.96 

.7469 

•1757 

124.^0 

.7496 

.1759 

124*23 

.8008 

.2274 

129.05 

.9520 

.6581 

152.17 

.9601 

.7058 

154.28 

.9844 

.8675 

161 .79 
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shall now calculate the activity coefficients and 
also estimate the inherent error in their values resulting 
from the use of experimental data. 


For calculation of the activity; coefficients, the vapour 
phase in equilihrium with the non-ideal liquid phase has been 
assumed to be an ideal mixture of non ideal gases. The relevant 

equations for the activity coefficients will then be 

(B^ ^ -Y°'^) (P-P^) + (2B^ 1 -^22^4^1 




P°x^ 


exp 


ET 


and 


PVo 

T, = 

^ -nO 


exp 


B22-V°^)(P-P^) + (2B^2-®11-®22^y?f 


K 


j 


(1) 


( 2 ) 


Equations (1) and (2) require- , apart from T-x-y data 
( i) the isobaric pressure P, (ii) pure component vapour pressures 
which are estimated by the Martin equation for methylcyclohexanol 
and by Antoine equation for methylcyclohexane , (iii)‘ second 
virial coefficients B-^ and B 22 which havco been predicted by 
emplo'i^'ing Kreglewski’s double square well potential and ( iv) 
the pure liquid molar volumes 7°^ and calculated from 

the analytical expression given by Yen and Woods. This 
equation requires the values of critical compressibility, critiea 
temperature T^ and critical volume V^, These critical constants 
have been obtained from literature for methylcyclohexane and 
predicted for methylcyclohexanol. 

The expression 23-1 2 -B-I Y-B 22 iias been assumed to be 
negligible compared to the other terms indicating that the 
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interaction virial coefficient is considered as an arithmetic 
mean of the pure component virial coefficients. This is the 
situation for an ideal mixture of nonideal gases where the equation 
of state for the nonideal gases can be given by the virial equation. 
Details regarding the use of pertinent equations for calculating 
t’'’e necessary information are being given below. 

(i) Antoinne equation^ for mcthylcyclohexane 


log-iQ P = A - B/(t+C) (3) 

Constants A,B and C have been obtained from literature. 

1 2 

(it) Martin equation for vapour pressure of methylcyclohexanol 


log-] Q P = A +■ B/T + CxT + DT^ (4-) 

+ ET^ + E log^Q(CT-T) 

The constants have been obtained by fitting the experimental data 

collected by A. ^dAklierjee of our laboratory. 

10 

(iii) Kreglewski jticdel gives the following equation for virial 


coefficients and B 22 ' 


B 

FP 


= 1 _ (p^_l)(e VkT_.,^ _ su/KT 


Y-'here R-| and R 2 S-r® given by the following expressions; 

1 3 


= 1 1 + 


Hr 


I 


s ■! 

Dr)T73j 


1 +( — 1 


V3'l 5 


where 7* is the liqiiid molar volume of methane at reduced 
temperature = 0.6| S is a constant = ,21 rounded off for 


-1 ) 

(5) 

(6) I 

i 

all substanc-,-4 

I 
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and Y* is the liquid molar volume of the compound at = ^6, 
The potential energy term u/kT is estimated by 


IL. - P / V*n1/3 

“ kT “ J c T '' V* ^ 

0 

is the critical molar density, T^ 


( 8 ) 


the critical temperature 


and T is the temperature at ■which we wish to calculate the 


virial coefficient. The packing factor b is estimated by 


b = .77 + ,1325 ( 9 ) 

The required values arc J"’q, T^ and Y* . Por methylcyclo- 

th 

hexane these have been obtained from the literature and for 

me thylcyclohexano 1 they have been estimated, 

1 1 

(iii) Yen and Woods equation for liquid molar volume 

= 4r^ = 1 + A (1-Tr)^'^^ + B(1-T^)^/^+h(1 

^ ° ( 10 ) 

is the reduced molar density 

T^ is the reduced temperature 

f _ is the critical molar density 

A,B and D are constants which are functions of critical 
compressibility alone. 

Then 

Y°^ the liquid molar volume is = 1 /(j’j.J^q) 

The critical constants for methylcyclohexanol have been 
estimated as follows; 

(a) G rit ical Tempera t ure; 

Tq = T-j^ (normal boiling point)/© 


( 12 ) 
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= .052 +.024 + *082 + *02 
= .178 

e = .567 + .178 - (.178)^ = .713 

Tq = 441 /. 71 3 = 620°K 

(1^) Or it leal Volp-tne ; 

Schuster’s method has been adopted using the lydersen’s 
volume increments as tabulated in Table 2.12 of Sherwood and Reid. 
Errors of 2 to 5"^ are expected. 

= Base value + (CH^)- 
+ 4(-CHp-) + 2(-CH) + 0H‘ - 

I 

in ring in ring in ring 
= 40 + 55 + (4 X 44*5) + (2 x 46)- = 383 cc/g.mole 
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C critical P ress ure ; 


Lyd e r se n * s me tho d g ive s 

■P = M 

° c# + 


(H) 


where M = molecular weight = 114 for CyH^-^O and is additive 
constant which is determined from molecular structure contri- 


butions as given in Table 2.6 of Sherwood and Reid. 

f 

1 = (CH^-) + 2 (-CH) + 4(-CH2-) + (0H-) 

= ,227 + (2 X . 192 ) + (4 X .184) + .06 
= 1.407 

= — Pii— , = 37.4 atm. 

(1.407+.34) 

Critical Compressibility 

7 — c c 

° = RT, 


(15) 


82.07" 


.281 


Estimated virial coefficients of Methylcyolohcxane are 

given in Table JEI and Methylcyclohexanol given in Table III,. 

To obtain an analytical expression for the variation of the 

virial coefficie nt with temperature a polynomial of the form 
2 

a + bT cT + ... has been fitted and the degree of polynomial 
is decided by subjecting it to E-test with 955^ confidence limit. 
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TABLE III 


VIRIAL COEPFICIEILPS OP MET I!TLCYCLOHB JAKE AKD METIJYCYCLOHEXAKOL 


fusing 

-u/ET-fc^ '^-TiTr 



L =.77+.^ 32, 5; (51^) 

in Kreglewski model) 

in°C 


^7^14° 

575*16' 

-1501 .56 

-5787.34 

578.16 

-1446.27 

-5550.97 

583.16 

-1394*19 

-5329*8n 

588.16 

-1545*07 

-5122,50 

595*16 

-1298.67 

-4927*90 

598.16 

-1254*78 

-4744*94 

405*16 

-1213.22 

-4572 .69 

408.16 

-1173*81 

-4410.28 

415.16 

-1156.39 

-4256 .94 

418,16 

-1100.83 

. -4111 *98 

425*16 

-1067.00 

-3974.77 

426.16 

-1054.77 

-5844*74 

453*16 

-1004.04 

-3721*36 

458.16 

-974*71 

-3604. 18 

433*16 

-946 *70 

-3492.75 

448*16 

-919.90 

-5386 ,6 c 

453.16 

-894.26 

-3285 .61 

458.16 

-869.70 

-3189*22 
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TABTiE 17 


EELEYMT INFORMATION FOR GAICULATION OF 
ACTIVITY COEFFICIENTS 


Antoinne Constants for Methylcycloliexane ; 

A = 6.82689 
B = 1272.864 

C = 221 .630 

Marfin Constants for Metbylcycloliexano Is 
A = 1 .6235139 
B =-943.73925 
C =-.140738 X 10"^ 

D = .3470786 x 10"^ 

E =-.41477109 X lO”"^ 

E =-.9089265 
Ct =627.98999 

Critical Constants for Metbylcyclohexane 

Z, = .251 
I = 572.1 “K 

V = 344 cc/g.mole 

V 

Critical Constants (estimated) for MetliyLcyclohexanol; 
= .281 
= 620°K 

V„ = 383 cc/g.mole 
P = 37.4 atm. ' 



(Table IV continued) 
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Coefficients in the polynomial fit for 'virial coefficient of 
CyHi ^ which has a, 5 degree polynomial ' . 

A = - .93186538 x 10^ 

B = + .91577548 X 10^ 

C = - .37652485 x lo'^ 

D = 4- .79423174 X 10 "^ 

E = - .8521 4885 x 10“^ 

E = + .37029944 X 10~^ 

Coefficient in the polynomial fit for virial coefficient of Cr^H^ ^^0 
which has 3 degree polynomial: 

A = -.98109325 X 10^ 

B = +.55025579 X 1-0' 

C = -.10869500 X 10^ 

D = +.73802514 x 10"^ 

Expressions for Constants in Yen and Woods Equation: 

A = 17.4425 - 214.578 + 989.629 - 1 522.06 

B = -3.28257 + 13.6377 Z^ + 1 07.4844 Z^ - 384.2111 Z^ 
when Z^ .26 

B = 60.2091 - 402.063 Zq + 501.0 Z^ + 64 I.O Z^ 
when Zq .26 


D = 0.93 - B 
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Error Analysis s 

As has been indicated earlier that rapour composition 
when regressed, on liquid composition gives ^ standard deviation 
of .006 v/hicli is greater than ,0037 the accuracy of measurement, 
and the equilibrium temperature T vjhen regressed on x gives a 
standard deviation of ,22°C which is again greater than .1°C 
the accuracy of measurement by copper-constantan thermocouple < 
Because of both these factors we keep the error in x as the 
fundamental error and the- errors in y and T are dependent errors 
and are calculated as follows; 


A rn dT 


( 17 ) 


^ V d. 0? 

The derivatives ^ and ^ are obtained by differentiating the 
respective polynomial expressions for y and T. To give an idea 
of the e,rrors in y and T these derivatives are being given in 
Table . It is clear from the table that because of steep ascent 
in the pure region of methylcyclohexanol errors incurred will be 
more. The errors so incurred will affect the calculation of the 
activity coefficients. It is our aim now to analyze the method 
of error propagation and see how ultimately t.hoy are reflected 
in In ? la "^2 

According to the general formula of errors if we have a 
quantity IT given by 

1 = f (u-j j U 2 . . . . ,Uj^) (18) 

then ^ U the absolute error in N is given by 
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TABLE V 


DERIViTIVES OE VAPOLR CQLIPOSITIORS AITD . 


EQTJTRIBE.ILt>l TET'IPERATIJRE WITH 

RESPECT TO 

X 

LIQ,UID COIffOSITION 

dx 

dx 

0.0505 

0 .0606 

8.1 41 1 

0.1234 

0,1309 

7.2224 

0.1882 

0,1865 

12.81 55 

0.2556 

0.2443 

23.7648 

0.-3177 

0.199'! 

28.3267 

0.5417 

0.2531 

33.2051 

0.5471 

0.2639 

33.4912 

0 . 64I 6 

0.4415 

55.1234 

0.6443 

0.4468' 

55.8951 

0.7469 

0.8069 

84.6485 

0.7496 

0.8229 

85.7489 

0.8008 

1 .2148 

99.7107 

0.9520 

5.4226 

■ 251 ,1919 

0.9 601 

5.8923 

271 .4175 

0.9844 

7.5809 

346.6897 
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As Au 

And since these e 
take the absolute 


+ 


'dUr 


A U2+ 


■ A 


A 


u 


k 


k 


rrors can act in opposite directions one must 
of the errors contributed by each, independent 


variable . 



( 19 ) 

( 20 ) 


A ( In 1 ) 



( 21 ) 



neglecting the error in the non-ideal term exp u^ . Ay is 
obtained from equation (16). /liP is the error in measurement 
of total pressure which directly depends on the error incurred 
in measurement of ebullio metric temperature of the reference 
liquid benzene. Corresponding to +.1°C at normal boiling point 
of benzene = + 2.1 4 mm. 



dP° 
= ^ 


At 


Since P° corresponds to methylcyclohoxano 1 the derivative is 

obtained by differentiating Martin aquation with respect to 

temperature. At is obtained from equation (1?). 

A similar analysis for error in A(ln'^p) is adopted 

dPg 

except that now is obtained by differentiating the Antoinne 

equation. 
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Y”-] _■ 

Since In = In * -i - Ini p 

^ 2 j. 

1 1 

absolute error in In -rq— is given by 

'2 


IS (In 


li' 

Tz' 


AdnT^)) + l^dn 


The values of in , In ‘'r2 s>^-^ ln'^-|/'’f2 
from exporimental data are reported in Table VI, 
errors incurred in each. See Pigure 8 . 


calculated 
along with 













CHM’TER 5 


THEBIvIODYUMIG CQlISISgEHCY TEST 

The G- ibbs-DuEem equation interrelates activity coefficien"^ 
of all components in mixture. Tberefore, activity coefficients 
calculated from thermodynamically consistent experimental data 
should obey the Gibbs-Duhem equation. If thev do not then 
it is concluded that the data are not thermodynamically consistent. 
This logic is the basis of all thermodynamic consistency tests-. 

The Gibbs Euhem equation for the binary isobaric data 
Is given as 

^ '^1 ^ "^2 ^ ^ n n ^ 

’ dx-j ^2 dxi dx-] 

This indicates that in addition to activity coefficient data 
one needs the excess heat of mixing data also. This is. a 
serious limitation since the excess heat of mixing data are 
rarely available, but it is a limitation which is unavoidable 
and pertains to all tests for thermodynamic consistency. There- 
fore, at best the test can only be approximate, Herrington 
has given a semi-empirical test for isobaric data and the heart 
of it lies in an estimation of the last term in the above 
enuation. 

Integral or area test of Eedlich-Kister and Herrington; 

This approach tests the data as a whole rather than 
iDoint by point and it has the important advantage that it tan 
be easily carried out quantitatively. 
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E 

^ - 
RT ” 


x-j In ? -j + X2 In f2 


( 2 ) 


Differentiating equation (2) 
d(g®/RT) ^ ^^*^1 

u tr^ -EpAn -. . /i ll .. — •y - I- - - I , t y 

dx^ ^1 dx>[ ^ ^2 


d InT. 


dXi 


+ 


In 


ln~T 


= In 


"fi 


1 


1 


Ji® ^ 
RT^ dx. 


1 


(3) 


On integration of equation (5) 


1 


=J in— '-ax 


E 

T 

“ 3 x 


1 


"1 


o 


I' 

■^2 


h' 


E 


1 


,.J 

o 


RT 


dT 


E 

S™ 

RT 


_E 


'1 


=1 


RT 


= 0 
x-j =0 


(4) 


Since excess free energies for pure components are 
necessarily zero. 

Therefore if one plots In -rrr- vs xa and neglects the 

I2 

heat of mixing term then the area under the curve should be zero. 

Eor the experimental data obtained for the system methylcyclcj 
hexanol - methylcyclohexane the area is negative, (Eigure Q). 

It makes no attempt to equalize to zero. In such a case even 
applying Herrington's test has no validity because the factor D 
which is defined as 

^ a rea above x-axis - area belov'/ x-axis 

^ arVa 'above x-axrs + area below x-axis 

turns out to be negative, equal to -IOO. 

E 

Actually for hydrocarbon-alcohol systems the g function 
behaves irregularly like for Benzene-BuOH system it behave's as 
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shown be low 5 


E 

One can then imagine the behaviour of h which has more assymetr]?’ 
E E 

than g . It is likely that h for methylcyclohexanol - methyl- 
cyclohexane system is highly assymetrical and irregular. There- 
fore in the absence of h data our test for thermodynamic 
consistency is not valid. 




THEOBETICAL MALTS IS 


In our experiment we Lave collected isobaric data for 
the system methylcycloLexanol-methylcyclohexane at 760 mm Hg. 

How does this data lend itself to an insight into the behavior 
of liquid mixtures? At first glance it does not offer us much 
hope for isobaric data is directly useful in design of distilla- 
tion units and it is the isothermal data which offers a chance 
for direct study on the behaviour of liquid mixtures. This 
isothermal data is more difficult to collect experimentally. Can 
we then predict the isothermal data from the isobaric data? 

This is a weighty question with profound implications. 


If we can, then we con predict P-x-y at different temperatures 
thus gathering information for all P-T-x-y. Let us form an 
image of what this means- 
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The answer to such an interesting question is in the 
positive, and it covers zeotropes as well as azeotropes. 

The structure of our prediction of isothermal vapour 
liquid equilibrium is built on the Eedlich-Kister expression 
for excess free energy, which is based on the conformal solution 
theory, 

E 2 , 

g = x^X2sl + BCx-j-Xg) + C(x-|-X2) +»..*■ (1) 

E 

Yt^ith B and C equal to zero, the expression for g is the same 

as given by regular solution theory, Bor strictly regular solution; 

the constant A is independent of temperature and S^, the excess 

E 

entropy of mixing, is zero and h -the excess heat of mixing, equal 

to . However the available experimental data shows that is 

E E 

not necessarily zero end h has more assymetry than g .. This 
indicates that A is not independent of temperature. Same reasoning 
would then apply for the higher constants in the Redlich-Kister 
expression. 

It was Haase, who first found the temperature dependence 

of constants A and B for cyclohexane-Methano 1 system, The expert- 

E 

mental curve of g for this system is represented by the equation 

g® = XiX 2 (A + Bx|) (2) 

in which A and B are the following functions of the temperature 
A = 2070 - 1 .85 T (3) 

B = -3107 + TOT (4) 

Haase then used these equations to compute S^ and h . However, 
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we can use the above explicit forms of A and B for our purpose 
of predicting the isothermal data. Calculate A and B at different 
temperatures and plug the obtained values in equation (2), Prom 
the expression of g one can obtain the activity coefficients 
vjhich can be used to predict P-x and y-x curves at the choosen 
temperature. The problem then essentially boils down to the 
determination of temperature dependence of the constants. This 
procedure is intricate and we shall develop the theory bit by bit. 

Malesinski has employed the regular solution theory l.e'. 

E 

g = Ax-| ^Xg and assumed A as independent of temperature. The 
value of the constant A has been determined from one experimental 
T-x data at x = .5. , This value of A has then been used to predict 
isobaric data over the whole rang^J of concentration.. Satisfactory 
results have been achieved by him for systems Like methanol-acetone, 
chloroform-acetone 5 heptane-beznene . 

However, the value of A can be improved by a process of 
successive approximations. If we take Tejjp-^pss a function of 
A and expand it by Taylor's series about the point then 


T = 'T' 

exptl “calc 


+ 



A + 




• • » « 


or ^exptl“^calc ^ fl truncating the higher deri- 

vat ive s 

The aim is to reduce the error in ^ T by the least squares 

C 

technique. Application of the least squares method will give 
the incremental value in A. A nov/ alters to A + ^ A, With this 
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new value of A, T _ over the whole range is again obtained. 

cl C 

In this manner, one can iteratively improve reduce 

the error in ,1 i.e. (T ^ -T , , ) . However a stage will be 
reached when no improvement is possible, and yet deviates 

largely from a^re still essentially dealing with 

■p 

one constant for expressing g for a strictly regular solution- 

E 

which is a symmetric curve for g when plotted as a function of 
composition. 

This is the stage when we must incorporate more number 

E 

of constants in the expression for g . We can do this in two 
different ways - (i) make B & C in equation (1) non-zero retaining 
their temperature independence or ( ii) make A dependent on 
temperature, this dependence may be linear or quadratic to fit the. 
experimental data, and keep B and C zero. Once these temperature 
dependent constants are obtained we can employ them for isothermal 
prediction. 

Let us novi/ proceed to obtain T-x- and y-x. 

The equations for the isobaric vapour-liquid equilibrium 
when a real liquid solution is in equilibrium with an ideal v.apour 


are 




= '^1 "^1 


(5) 

Py 2 = '^2 ^ ^2 

The equation for the boiling temperature isobar is obtained by 
adding (5) and (6) and dividing by the total pressuire P. 

l.e. =1 (7) 
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where \ = Pressure of 1 at Temperature T 

i Total Pressure ' 

The equation for the vapor composition in terms of liquid compo- 
sition will then he 


= 


r ^ r 


( 8 ) 


'^i 

The activity coefficients appearing in equations (7) and (8) are 
obtained from tiio following relation 


f'. = L M? = 1 


^(nrpg®) 


In - ^ ^ -J-W: 


’ ^T “'2?i 


(9) 


Equation (7) is implicit in T and it does not lend itself to 
easy calculation of T, This difficulty can be overcome by using 
another form of the equation for the boiling temperature isobar, 
which can be derived by using the excess free energy of mixing. 
Moreover at this point we wish to incorporate the non-ideality in 
the vapour phase also. So that 

2 


Pv 


1 


exp 


B^^)(P-P^) -5^2 ^2 ^ 


RT 


and 

or 

and 

where 


j . (V^^- B22)(Mp -<^12 

'^2 ^ RT 


Py2 = P22:2 ? 2 
Py^ = P® x-| i 1 exp u^ 

Pyp = P2 ^ 2*^2 ^2 

{Yf^- B.,|)(P.-Pp -^^^2 


( 10 ) 

(11) 

( 12 ) 

(13)- 


u 


1 


RT 


(^2^- ®22)"*^^--2^ "Sjg yp 


V.r 


RT 


i.CEH 


A .46|'p 


^ , 
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V" 


C'-'j2 “ 2 “ 1 “ ®22* 

The excess free energy of mixing is given by 


E 


RT (x-i In 


j_ii I -| + ^2 In Y 2) 
Using equations (12) and (I 3 ) 


(U) 


E r -^1 

g = RT (x. In 


P° x^ exp 


+ Xo In 


Py. 


u 


1 


^2^2 ''^2 


) 


\ p p 

= RT i In — r + Xo In -rr- + In 

Pi Pg 

^2 


Pi 


+ 


X2 


In 


Xi Ui 


^2 ""2^ I 


or 


E 


RT ( V + w - u ) 


( 15 ) 

(16) 


where v, w and u are defined as 


P P 

V = x-i In — — + Xq In — r 

^ P? ^ Po 

y- 

w = x-i In 4 - Xo In —• 

I X-^ ^2 


la = x*| u-] 


+ Xo n 


2 ^^ 2 


(I7a) 

(17b) 
(1 7c) 


Prom eouation (16) 

^E ^ 

E(v + w - u) 

E 

' If the theoretical . expression for g is fixed the above 

equation can be used to obtain T. Equation (18) though implicit 
in T like equation (7) is not as unwieldy and T can be improved 
iteratively with the first approximation on the right hand side 
taken as T = T^Xi + TgXg. 
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However equation (18) belies the trust put in it. 

Computer calculations revealed that v is very sensitive to 
temperature and even a little off approximation for T on the 
right hand side makes the expression (v + w - u) negative , ''-ith 
the result that one obtains negative temperatures. To counter- 
act this problem and make the iterations converge we had to 
express the term v differently. 

Prom Claussius-Clapeyron equation, we have 

RT In - T) (19) 

^i 

if we assume the molar heat of vaporization to be constant over 
the temperature interval in question. Here is the pure compo- 
nent boiling point at pressure P and is the entropy of 

vaporization. Plugging this relation in equation (16) one obtains 

g® = x^ /\S^(T^-T) + Xg AS 2 (T 2 -T) + RTw - RTu (20) 
Transforming the above relation 


T _ T X -^i + T X ^ EIw Mu 

I +13X2 ^3^ -^^ 3 ^ -_^ 3 ^ 

where As^= x^ + X2 A ^2 


It is this equation which w« are going to use for getting the T-x 


relation. The first approximation for T is again T^x-| + T 2 X 2 , 

E 

The entronies are calculated from equation (I 9 ), g is given 
by Redlich-Kister expression, u is the expression for non-ideality 


of the vapor phase. It is the term w which has to be made more 


explicit. 
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yi T2 

w = In ^ + Xg in — 


c. ^2 

Tn terms of relative volatility ^ 21 

1 


^ ^ 

^1 ^ ^2 < 21 


72 


X 


2 ^^21 


X , 


2 X -^ + 2^ < 21 


w = X -1 


+ 


““s r " 

ln\ 1 + X2 ( <21 " ” ^2 I ^ + ^( <21' 

L, 


^2 In <21 

= X2 In <21 - In j ^1 + X2(<21 " M 

where the relative volatility {<21 defined as 

/ _n 


fV -^2/^2 

^■21 yi /xf 


Pg r2 exp u, 
P° r,j exp u-j 


( 22 ) 


(25) 




•V *’^2 

and ln< 2 i= In “ + In vv- + U 2 - n 


•1 


tl 


(T-Tg) - ^ 


w as a function of In <21 alone is given by 


(24) 


w = ■- x^XgCln < 2 < ^1^2 

X -] Xg 


(ln<n^ 






(1 -6x-iX2)(ln<2l + X^X2(X^-X2) 


X -1 X 


(ln.£)(,2i)^ + -7^ ~ 50x1X2) 

+ 120x^X 2) (ln'rr<,2l)^ 
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Having obtained the T-x relation, we can also obtain y-x from 
equation (8) , 


where 




P.- exp u. 


1 P 

Having obtained a useful form of equation for T, which poses no 
convergence problems, we now proceed to our main objective of 
obtaining the constants in the Redlich-Kister equation. We shall 
consider the general case of quadratic temperature dependence of 
the first two constants 1 and B. 

,2 


A = + C 2 T + C^T 


B = C, + Cj- T + 

4 5 6 


(2 6a) 
(26b) 


so that 


E 

g 


— ( C-i + 


CgT + 


) x-^X2 


4- (C^+C^P + CgP^)x-jX2(x^-X2) . 

(27) 


Particular case of Maleisinski ' s predictions 

= C-| x^ X 2 (28) 

C-] is obtained from the experimental point at x = .5. Only- 
temperature value is needed, not the vapour composition. And 
since Malesinski 'has neglected the non- ideality of tho vapour 
phase the term u in equation (16) is neglected. Prom equation . (1 6) 
in conjunction with equation (28) 


C^= 4RP* (v* + w*) (29) 

*indicates that the values v and w are calculated at temperature P-^ 
corresponding to x = .5. Since Malesinski makes no attempt to 
improve the value of this constant C-| , no least square analysis is 
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.required. And therefore one essentially obtains predicted ' 

(T-x) curve employing equation (21 ) with u term equal to zero 
• « 
and predicted (y-x) from equation (8). 

More than one constant with temperature dependence 

included* 

C = (C^ Cg C3 Cg) ( 30 ) 

We will determine the elements in this linear vector C by least 
square analysis, and the values are improved iteratively. To 
find a first approximation we assume that the solution behaves 
like a regular solution, so that only is non-zero. C-| is 
evaluated at equimolar concentration for the working of the least 
square technique (T-x) experimental data sre required. 

4Ef* (v^- + - u*) ( 31 ) 

T is calculated for the experimental liquid compositions 
using equation (21). Pirst approximation for T = f x-] + 12 ^ 2 * 
Equation (21) is used iteratively?- at each liquid composition till 
no further improvement in value of T is possible. Then the 
temperature residual /.)> T = - '^exptl calculated for each 

data point. 

£\, T is a function of the elements in C. By a Taylor 
series expansion of '-^expil '^calc 

''exptl = ’^calc P" derlyatlvea. 

T=i ^ 

We shall neglect the higher derivatives and then 


A 


I 




^ 'Sg, 


M is the no. of elements in C. 

(33) 
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The first derivatives of T with respect to elements in C are 
ohtained as shown below. 

Let us say that 


T 


f (T, C^) 


i = 1 to M 


(34) 


Then 

and 


2 l 

?jT 


df 

^ - M 
ac^- dc. 


"7^ f 

dT +J^ dC, 


"df dT . 2 ) 


: T dC^ dC^ ’ dF7 " ^ 


df/dc. 


dT 

-H) 


for i (35) 


(36) 


where f is given by equation (21). 

to 

In equation (36) we have yet/analyse the partial derivatives 

9 


f 9 f 

— and We shall now present the analysis for these 


4)C 


1 


9 T- 


partial derivatives. 

ni_f 
'~o T 


Analysis for Taking the derivative of equation (21 ) ^ 


^ f 
■9 T 


r 

= T-jX-j -rf + T2X2 


d-T 


JLe 

T 


+ RT(^| - “ A 


9 


T 


d.As. 


dT 


d As° 

Jt 


Tt^ -t) 


d In p: 


R In 


RT 


1 


dT 


+ 


As, 


(37) 

(38a) 


a similar equation will be obtained for the second component. 


o 


d In P7 


Por 


dT 


v\?e employ usually the Antoine equation. 


A 


= X. 


A 4 

ai 


^ ^ "dT 


(581*) 



58 


From equation ( 27 ) 


'1 .E 

i-£_ 


= (C^ + 2C^ T) x-|X2 + (C^ + 20 gT) x^X2(x-|-X2) 


Ihe temperature derivatives of the term u are intricate. 

Im ^ 2i«2 

T 1 v T ^15 r 

-i 


(38c) 


{38d) 


^i= 


(vf - B.^)(p-p“) -d 12 yf 


ET 


In the above expression liquid molar volume, virial coefficient 
vapour pressure of pure component i, are all temperature 
dependent. We shall neglect the temperature dependence of 0-^2 
Prom the expressions of vapour compositions y-^ and 
in terms of the relative volatility ^21 have 


<3y2 

_ ^2 ^1 

^ C 2 I 

"dT " 

-0(..,21 

ai 

<3y-i 

2 ^2 

ao('2i 

"df " 

"^1 x^ 

dT 


(38e) 

(38f) 


Using equations (38e) and (38f) we get 




“1 




I 


u. dV?^ dB. 

+ ( -- - -^— ) (P-P^ )/ET 


■( 


~ ^ 

“ ''PI 


d:^ 

"dT 


(38g) 


2P O 


yg yi 


12 


21 


/m 


' d^2 


Vl 


2 . “22, ,o 


d 


+ (: 


I ‘ MT dT ' '^“"2 

(Yg^- B 22 ) (3^2 2pS^2 


RT 


dT 


+ 


) (P-?p)/RT 

3 

^1 X. 


^2 ‘^'^21 


1 


dT 


(38h) 
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From equation ( 24 ) 


21 

T' 


= <21 

• CvJ 

1 — 1 

! 

'3 T 

. I 

id In P? 

= ^2P 

I 41 

1 dT 

3 

Ug '3i 

^ T - - .T 


d Ln P 


dT 


1 ■ 3 gE 

PT // ^2 


RT 


Xj'' 


1 

rn 


(38i) 

Equations (38g), (38h) and (38i) can Toe solved simultaneously to 
get explicit relations for 


'3 U2 


O 

u 


1 


From equation (24) we get the expression for 


■9w 


'7)w 


r -/ (x-i-x^) ./ 2 

x^Xg 2ln'^^-^ + — 2 (FnC\,2i) 


+ 


(1-6 x-^ X2 ) 



2 2 

1 - 30x^X 2 + 120 x.^X2 


^ — ( ln''^■^^21 ) ' ■*■ 


^ (X- 1 -X 2 ) (1-1 2x^X2 




'\21 


^4 

)5 + ... 


Un ^ 21 ^ 


1 i^\21 

0(2^ 


d T 


(383). 


Analysis for 


Pf/'B 


Taking the derivative of f with respect to constant C. 


P 


1_ 'U&L * M 


_ CO _j_ ■ - - . 

"FR ° FP AS, 


w 


// 




RT O u 
iSo'?C. 


(39) 


3 


u 




^Ui 




0^2 




It is only the -R ■\2, term in expression for u which, is 
function of the constants C^. 


a 



6'0 


From the expression for' u-| , Vl 2 and equations (38e and 


(38f) 


('Juo 2p (A 


2) ^i 


p'"" 

12 2 
E.T ^1 ^2 


■ ln{y2'i 




r 
2P 6 


12 


JXn 


■'7 


C, 


RT ^2 


21 




(40a) 

(40b) 


From equation (24) 
3 In 0( 




'•_21 
i 


Jun 


d 


u 




i 9^- 


1 i__ 

^ RT 




■) 


(40c) 


1 — ry 2. cJ 2 

The last term on the right hand side in the above equation is 
determined from equation (27). 

Equation (403), (40b) and (40c) solved simultaneously to 


9 


u. 


6-'^1 


get explicit relat ions f or —• rrr” 5.— —7^ 

tj ^i ^l 

id ? / X 2 ) 
O 


and 


{) ^>^ 0(2 

-^“c- 


1 


0 In ::X2i 


- 2^522 El E2 


(40d) 


2 


E 


20, 

1 to M. 


- can easily be obtained from equation (24) for all i from 




The expression for ,.-4^ will be similar to .-X-^ as 

O^i -Q In 0^4^/ 

ffiven by equation (38j) except that instead of y-Tfi — ; 

u i 


/) 


\ln.:5^ 


.,21 


O ^i 


will be used . 

With the help of e quations (37) to (40) we can no?; 
determine the first derivatives of temperature with respect to 
elements C, in C. Returning back to equation (33) 


4^ = Z ^ ^ 


33) 
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■''e determine the changes in which will most nearly 

reduce the temperature residuals to zero. Least square technique 





M 






1 

M 



d T ?) T 

Cl, 


n 

a M 


+ 


b c 


2Z 


('J T 3 


a T 
0 


1 


2^ 


+ . 


G 


T^2 

M 



' 3 T 



A simultaneous solution of the above set of equations (41) yields 

,'■'5 C^, Adding these increments to the initial values of constants 

v/e find improved approximations. Since the second derivatives 

of T with respect to are not zero, these are not necessarily 

the best values, and it mav be necessary to repeat the process 

until do not change significantly. This procedure has not 

utilized the (y-x) experimental data, (y-x). relation can be 

predicted for the isobaric (T-x) data under consideration. 

The above process essentially gives us the constants A, 

E 

B, etc. in the Redlich-Kister expression for g , as functions 

of temperature. Erom the values of C-, to 0^^ the constants 1,B 

E 

can be reconstituted at different temperatures, g is then 
obtained as a function of liquid composition alone. One can then 



CHAPTER ■ 7 


DISCUSSION OH THEORETICAL KESIJDTS 


In order to e valuate the method developed in Last 
chapter (Theoretical Analysis) for the prediction of isothermal 
data from isobaric data, predicted results have been compared 
with available experimental data. The comparison has been done 
both in tabular form (Tables VII and VIII) as well as by 
plots (Rigs. 11 to 27 ). 

The systems chosen for comparison are non-polar + non- 
polar mixtures, non-polar + polar and polar + polar mixtures.. 
These belong to the azeotropic as well as azeotropic groupings. 

As a measure of deviation between experimental isothermal 
data and predicted isothermal data two quantities have been 
defined. 

( i) Stdp, the standard deviation in total pressure 


n^ 

2 ^- ^^calc"'^exptl^ 



and ■ 

( ii)Stdy, 


the standard deviation in vapour phase composition 
n 

^^calc ^exptl'^ 

* i 

Stdy = ~ — 


n -1 


where n is the number of data points. 

The extent of non- ideality of the system is given by the numerical 
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'TABLE TII 


PREDICTED 

ISOTHEBIAL DATA AT A 

GLANCE 


Temp • 

system “K 

Std. dev. 
in P 
in mm Hg. 

Std •dev. 
in Y 
in mole 
fra ct ion 

Numerical 
value of A 
in Cal/g.mole 

Form of the 
constant 

I. CCl^(l)+Cyclohexane(2) 283^^ 
(non-polar+non-polar) 31 3^® 

1.74 

2.06 

.0067 

.0032 

78.83 

78.83 

Independent of 
temperatare 

C (1)=2.238, £,(2)=2.220 

553^' 

3.69 

.0039 

78.85 

IT 

P = 760 

345^^ 

4.98 

1 .1 7 

.0032 

.0007 

78.83 

55.13 

ir 

A=-958.5+2.95T 

II.CCl^+CgHg 

(non-polar+non-polar) 

2.(1)=2.238, ’^(2)=2.284 

S1 

P = 760 mm. 

31 3^^ 
343^^ 

.51 

.96 

.0027 

.0017 

65.51 

, 65.5T 

Independent of 
Temperature 

II 

(pig. 11 and 12) 




' 


III. Benzene+Heptane 

( nonpolar +nonpolar) 

353^^ 

5.97 

.0057 

249.16 

II 

t(l)=2 . 2847^(2) = 

P = 760^^^ 






IV. CHClj + CgHg 

(nonpolar +nonpolar) 

36 

298^ 

14.65 

.0559 

-371.78 

11 

;4(1)=4-8065 4- (2)=2.284 
P = 760 mm 






V. Acetone+Benzene 

37 

318-^ ‘ 

9.72 

.0114 

189. 5 A= 

=1 309. 37-3. 5T 


g(l)=20.7, g(2)=2.284 

^ 38 
p = 760^ 

(Fig. 14) 
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(Table YII continued) 


TI. CCl^+Acetone 318 

£(l)=2.238,g,(2)=20.7 

39 

P = 760 mm 

(Fig- 15) 

YII . Ethy la cet^te-f Propanol 315' 

£(1)=6.02,£(2)=20.1 

U 7^0 41 

P = 760 mm 

(pigs. 20 & 21) 

YIII -Ethyl Ether +Ethanol 

e.(l)=4.555,P(2)=24.50 

42 

P = 760 mm 
(Figs. 17, 18 & 19) 


40 


6.22 


.0147 


495.71 


IX. Acetone+Ethanol 

p (1)=20.7,2(2)=24.50 

^ '45 

P = 760 mm 

(pigs . 15 & 16 ) 

Azeotropes 

X. CgEg +rX 3 gH ^2 

£,(1)=2.284,^(2)=2.220 345 
48 

p = 760 mm 
(Pigs. 11 & 22) 


51 5^' 

3.40 

.0215 

607.17 

555^^ 

6.59 

.01 34 

512.23 

27545 

5.50 

.0134 

709.73 

283^^ 

8^8 

.0151 

71 2 .42 

295 ^^ 

' 10. -72 

.0125 

71 5. 81 

298^^ 

9.40 

.' 01 76 

717.76 

303 "^^ 

15-89 

.0135 

719.88 

515^5 

19 . '22 

.0142 

724.64 

505 ^^ 

2.13 

.0130 

511.00 

515 ^^ 

1.66 

.0066 

437.66 

5.45 

.0109 

495 .82 

321 ^^ 

.90 

.0055 

437.66 

6,89 ; 

.04-13 

| 80 .65 

315 ^^ 

.90 

.0022 

272.60 

) 545^^ 

1 .68 

.0028 

272.60 


A =-5402.46+20.271 

-.025 T^ 

A =-2095-56-4.751 
A =-2095.56-4.751 

A = 903-48r1 .'651 

+.0035 1 

n 

!! 

!1 

W 

M 

A = 1-089. 74-1 .90T 
Independent of Temp. 

A = 1089 . 74 - 1.901 

Inde pendent of 1 emv 
A =1089.74-1 .901 

Independent of Tem] 
u . 
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(Table VII continued) 

XI. Methanol- 

Ethy la estate 

£.(1)=32.63, £(2)=6.02 
P = 760 

(Figs. 23,24 & 25) 

XII. Ethanol -Benzene 

£,(1)=24.30,^(2)=2.284 

53 

P = 750 inin'^ 

(Figs .26 & 27) 


513=“ 

6.56 

.0167. 

3255° 

3.06 

.01222 

333^°1 

7.73 

.0132 


3135^ 

9.37 

.0557 

328^^ 

11 .93 

.0395 

333^’ 

17.65 

.0348 


679.70 

643 .29+ .121 

664.70 

Independent of 
Temperature 

4=663.92 

Two constants 

B=-14.90 

Independent of 
Temperature 

1094.45 

1=1885.32 - 
2.53T 

1204.82 

1=1 25 62-. 97 - , 
39.O6T+.O74T' 

1099.87 

I! 


£(i). £(2) are the dielectric constants of the , components# 
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Ca^PARATIVE STllLY OF STAIiDAHD DEVIATIOIjl IE PRESSURE 


Standard deviation in pressure with 


System 

Single 

constant 

Linear temp, 
dependence of 

A 

Qmdratic 
temp.dep. 
of A 

Two constants 
temperature 
inde pendent 

CHGl5+CgHg(298°K) 
(nonpolar4- nonplar) 

14.64 

45 .65 

48.43 

21 .83 

CCl^+A c et one ( 5 1 8 ° E) 
(nonpolar + polar) 

12.17 

12.02 

6.22 

11.69 

Ethylacetate + 
Propanol (51 5°K) 

8.17 

3.40 

20.99 

7.18 

(nonpolar + polar) 





Acetone+Etbanol(505°K) 
(polar + polar) 

5.53 

2.13 

7.OT 


CgHg+Cyclohexane(31 5°K) .90 

(nonpolar+nonpolar 

azeotrope) 

9.05 

9.29 

4.25 

C 2 H^OH+CgHg( 528 °K) 

30.76 

18. 42 

11.93 

26.19 


( Polar+non-polar 
azeotrope) 


Exptl. data o— 

Prczdicted with a Ct 

s (fig le cohstan t;^^^ ^ 


at 70*' 




C6H6-CgH^2 at70°C 






— • Exp tl. data 

V Predicted with 
single constant 


i I Icju id compos i t ion 
‘^oedloride + Benzene. 










1or Accton(?(] )- Ber^rt^r.e 


Exptl. data /o 

Predic ted Wit h. iinetir c 

temperc 

iture dep- y/ P 

endece 

of Vv 

first c:or>stant'^/ ° cc 


■ y v i 


/ ^ 





/ 

V ■ 

/ 

oc„ , • 

/ XJ 


/ ^ 


/v 

o - Predic ted wi thou t 

o 

t e rr, p . d epend e n ce 
b\ the f irst constant 







Predicted with linear 
terrip. dependence of 
the first constant 


Exptl. data 




'hquid.ccmf: 
rtane r tt^.rjroL 







Predicted with 
quadratic temp, 
dependence of the 
first constant 
Exptl. data 


Fig. 17 -Vapor pressure vs. liquid composition for 
, -r/ (-sElfeyl Eth<?r-Ethanol. 




'Or composition vs. liquid composit 
ilthyl ^ther - Ethanol . 



o Predicted With 


quadratic temp. 


dependence o f 

^ 1 

r t,he first constant 








Exptl. data 

Predicted with 
linear temp, 
dependence of 
the first cons - 


T=40t 5 


Fig.21-Vapor>compositiori vs. hquid composition for 
, . Ethyl. Acetate - Propanol. 













Fig. 26 -f quitibrium curves tor Ethanol-Benzene at 40°C. 






Fig. 27- Equili brium data for Ethanol - Benzene 

at 
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value of the constant A. The extent of polarity ’of the .compound 
is gauged by the value of the dielectric constant. 

Before discussing the obtained results a basic contra- 
diction in the analysis will, be outlined. In order to accommodate 

the assymetry in isothermal g one must have more than one constant 

E 

in the Redlich-Kister expression for g , and in order to include 
the concept of irregular solutions (S ^ 0) these constants must 
be made temperature dependent. That is when isobaric data is 
analyzed one calculates more number of constants (0^^) than are 
apparently needed for isothermal data (A,B,C etc.). The increase 
in number of constants (C^^) will reduce the number of degrees of 
freedom and thereby lessen the chance of improving the values of 
the constants by successive iterations in the least square 
analysis. Therefore most of the results have been obtained with 
one constant A (implying regular solution) but with temperature 
dependence (i.e.-not strictly regular). 

looking at Table VII it is found that for non-polar 4-non- 
polar mixtures the behaviour is strictly regular- the constant A 
has no temperature dependence. The value of A evaluated from 
experimental isobaric data is valid over all temperatures. The 
predictions agree well for systems CCl^ + Cyclohexane( I) , CC1^+ 
CgHg(lI)j Benzene + Heptane ( III) , and Chloroform + Benzene (lY) , 

Bor systems (I), (II),' (HI) the standard deviations in pressure 
are less than 5 mm, and standard deviation in vapour composition 
less than 1^. That (IV) has much larger deviations in pressure 
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and vapour composition can ‘be due to its higher non-ideality 
(A = -37’1.78), slightly polar ^nature of chloroform. However 
from Table VIII it is seen that results for (IV) are best with 
one constant only, two constants (A,B) giving higher deviation 
in P. That there is a slight anomaly in CCI^ + cyclohexane at 
1 = 3430 K - better results are obtained with linear dependence 
on temperature of A, indicates that non-polar + non-polar systems 
can be slightly irregular. 

for the non-polar + polar systems - Acetone + Benzene (V), 

CGl^ + Acetone (VI), Ethylacetate + Propanol (VII) and Ethyl Ether + 

Ethanol (VIII) , we observe from Table VII that (i) nonideality 

is more (A is higher than in non-polar + non-polar mixtures) , 

( ii) there is a clear dependence on temperature indicating the 

irregularity of the system, (iii) the standard deviations in P 

and y based on single constant A are more than for non-polar + 

non-polar systems indicating that there exists more assymetry in 

isothermal g , however ( iv) irregularity of the system is a more 

E 

important consideration than the assymetry in isothermal g , for 
Table VIII indicates that for systems CCl^ + Acetone and Ethyl Abets,. 
+ Propanol stdp is less for a single temperature dependent ionstant 
than for two temperature independent constants. The «hanges securing 
in A values for _^T = 10°K at T = 300°K are 13?^, W° and. 6^ 

for systems V,VI, VII and VIII. That larger deviations are 
obtained for ethylether and ethanol with' temperature dependent 
constant, and there is negligible change in value of A with 
temperature indicates that perhaps for this system irregularity 
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concept is not that important, one could do better with higher 
constants temperature independent in Eedlich-Kister expression. 

In all these predictions one thing should be noted that these 
isotherms do not intersect the isobaric data, that is these 
predictions are essentially extrapolations. 

For the polar + polar system Acetone + Ethanol the value 
of A is temperature dependent but it appears that it is not linear, 
it might be exponentially decaying with temperature. Again tn this 
case the predictions are essentially extrapolations. 

For the azeotropes not many systems have been studied. 

For the close boiling, non-polar + non-polar mixture CgHg+CgH ^2 
the behaviour is regular, indicated by predictions with stdp <1.2 mm 
and stdy ,55^. Among the non-polar + polar systems - for methanol 
+ ethyl acetate dependence on temperature is .weak and one could do 
as well with two temperature independent constants; for ethanol 
benzene temperature dependent constant is better than just two 
constants, however the deviations in P and y are large indicating 
that two temperature dependent constants might have been better 
than either a single temperature dependent constant or two 
temperature independent constants. 

In this discussion it must be remarked that all the mixtures 

studied have had had a unimodal plot of g (T,x) vs. x. But there 

few systems which are not unimodal and have more than one peak - 

24 

for example Acetone + BuOH has one minima and two maxima and 
2 5 

CCl^ + CH^OH has two minima and one maxima. In such cases one 
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should '‘have a higher degree of polynomial expression for g^. 

For the systems Acetone + BuOH and CCl. + CH^OH it should be 

4 5 

4th order - having four constants A, B,C,D in Bedlich-Kister 
expression. 

For the system Methylcyclohexanol + Methylcyolohexane 

(CyH^^O + In " 1^/12 has three peaks - two maxima and 

/ ■ 

one minima, that is In ^ order polynomial in x. 

And since Ini-^/fg = VeT x ( 2)g®/3x-|), g® is necessarily 
fifth order. Hence, Malesinski's prediction with single constant 
will not suffice. 

COl^fCLlTSIOH s 

Prom the above results it is concluded that 

(i) non-polar + uon-polar mixtures behave as strictly regular solutions i.e. 

El E 

S =0, and therefore in g = 1 x^x^, A can be used for all temperatures. 

(ii) ncn-polar + polar mixtures show a clear dependence on temperature of 
the constant A, indicating that the nature of association alters with 
temperature. Same is the case for polar + polar mixtures. However, 
the dependence on temperature need not be linear or quadratic, it can 
have other forms of dependence - say exponential. 

(iii) Isothermal data can be predicted even at those temperatures which do not 

lie in the range of isobaric T-x-yi that is extrapolations are possible. 

E 

(iY)the shape of isobaric g gives a clue to the nimber of constauts 

E 

needed in the Redlich-Kister expression. For a g with n peaks, n 
constants are needed. 



CBAFTER 8 


RECaMEIffiiATIOIS 

For the system methylcyclohexanol + methylcyclohexane the thermo- 
dynamic consistency test remained invalid because of the absence of heat of 
mixing data. It is recommended that the heat of mixing data should be 
experimentally measured for the above system. 

In all the isothermal predictions carried out (i) experimental isobari 
data was used (ii) interaction virial coefficient was assumed to be an 

arithmetic mean of pure component virial coefficients ^ and B 22 single 

precision was used in Fortran IT on IBM 7044* It is recommended that the 
isobaric data should be smoothened and tested for thermodynamic consistency 
before employing for isothermal prediction, the interaction virial coefficient 
should be estimated from experimental data or predicted using a suitable model 
with realistic combining rules. Kreglewski’s model can be utilized for 
calculation of B ^2 ^ function of temperature, u ^2 being given as a harmonic 

mean of u^^ and U 22 3-1^4 packing factor b calculated as a function of pseudo- 
critical temperature. Double precision should be used in order to match the 
sensitivity of the variables. Calculations should be performed with more 
number of degrees of freedom, i.e. more isobaric experimental data points. 

To get better isothermal predictions the isobaric plane should be selected 
such that extrapolation is avoided. 
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APPEICDII 1 


details op THEMOCOTJPXE-C iOLlBRAIIOE 


" I sought a method which above all would be rapid 

and simple, and decided on the use of thermo-electric couples..." 

Henri Le Chatelier(1 885) 


Principle of Thermocouple 



This is a very simple thermocouple circuit employed in 
the measurement of temperature. According to the laws of thermo- 
electricity small electromotive forces are generated in such a 
circuit. Briefly these laws amount to the following information. 

(1 ) If each section of wire in the circuit is homogeneous, that 

is if there is no change in composition or physical properties 
along its length, then the emf in the circuit depends only 
on the nature of the metals and the temperatures of the 
junctions. 

(2) If both of the junctions involving a particular homogeneous 
metal are at the same temperature, this metal makes no net 
contribution to the cmf. 
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(3) If all junctions of the circuit except one are held at constant 
temperature, the emf in the circuit will be a function of the 
temperature of the remaining junction, and can be used to 
measure that temperature. It is customary to prepare tables 
giving this emf as a function of temperature for the case where 
the reference junction (or junctions) are held at 0°C.”^ 

The thermocouple is a fairly accurate and stable temperature 
measuring device with rapid response and low maintenance cost. The 
commonest thermocouples are chrome l/alumel, iron/consantan, cooper 
constantan, and platinum/platinum- rhodium. Selection of the type 
of thermocouple depends on the upper and lower temperature limits 
the type of atmoshphere surrounding the measuring junction. & 
present study copper/constantan thermocouple has been selected which 
is suitable for -200 to 500°C range and the approximate change in 
emf per 100°C is 4.3 m7. It may be mentioned here that copper/ 
constantan thermocouples are subject to continuous oxidation since 
the metals absorb their OYm oxides. For this reason such a thermo- 
couple should not be used in strongly reducing or strongly oxidizing 
atmospheres. In addition to self -contamination, a similar and 
much more rapid fallof emf will follow any exposure to metal vapors. 
Mechanical failure of the hot junction can occur due to over-heating 
or excessive heating. Copper has a positive thermoelectric power 
with respect to Platinum and Constantan has a negative one. Thermo- 
couples can be made by employing the use of two dissimilar wires, 
but manufacturing methods have now been standardized and normally 
one uses these manufactured tnermocouple wires. The manufacturers 
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follow the simple code of putting postive wire in a red sheath and 
the negative wire in a blue sheath. Sheaths are used to prevent 
any short circuiting between the two v^ires. Each thermocouple type 
has a definite emf temperature relationship and this need not be 
linear. Since thermodynamic studies involve accurate temperature 
measurement it is intended to calibrate the thermocouple against 
a standard Platinum Resistance Thermometer. This over-all cali- 
bration v^ill also take into account the uncertainties in the circuit. 

The uncertainties arising in the over all circuit will 
briefly be mentioned. Each measuring junction has its own peculiar 
characteristics depending on such factors as materials, methods of 
joining, age, and history of junction. Yet these measuring un- 
certainties are small and are accounted f or in overall circuit cali- 
bration. Reference (or cold junction) can introduce uncertainties 
of + 1°E. These can be minimized by using the same reference junction 
in the field as was used in the calibration. The uncertainties can 
be reduced by controlling the reference junction. In this study 
cold junction (copper-copper and constantan-copper ) has been kept 
at ice point. In order to reduce the errors that may occur because 
of the presence of air bubbles in ice and also to ensure exact melting! 
point, ice and water are used in a slurry mixture. Copper leads 
of gage 20, enemel coated are used to contect the cold junction 
of the thermocouple to the potentiometer. These wires cause no 
uncertainty when they are used between the reference junctions 
and the potentiometer. The potentiometer introduces no uncertainty 
that is not accounted for in the overall calibration. 
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It has been mentioned earlier that the calibration is to 
be done against a platinum resistance ■ thartno meter. Platinum is 
chosen as the material for use in resistance th'ermometers because 
of its stability, ductility, resistance to oxidation, availability 
in a highly pur© state, and absence of irregularities in its 
temperature vs resistivity curve. The sketch of platinuai resistance 
thermometer is shoym in figure. 



1 ! mica cross 5 cm. long and .5 cm. diameter 
B s mica spacers 
C % base with binding pasts. 

The actual thermometer with its leads is inserted in a glass tube 
and hermetically sealed to a base with leads. Compensation for 
the change of resistance of the leads v;ith temperature is most 
simply effected by installing a loop of tbe same wire as the 
leads between the base and the point w'here the leads are welded 
to the platinum wire. 

According to lltueller the essential elements 

required in precise resistance thermometry (are) a resister properly 
mounted and protected, a means, for measuring its resistance, and 
a relation between resistance and temperature The mounting 
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and protection has been described _ already. The means for measuring 
resistance that is employed is a resistance recorder chart which 
measures resistance upto throe decimal places.- third place being 
determined by eye estimation. The recorder employs the Wheatstone 
principle. As for resistance - temperature relationship, Oallendar's 
equation describing the variation in Platinum Resistance with 
temperature is still pertinent. The Callendar formula is as follows 

e*qq-b°' + ^ (w - e tto 

(A 0 

An alternate form which is completely equivalent to the above 
expression is 

(1 + At' + Bt'^) 

A and B are related directly to and 
A = ( 1 + ^/lOO) 

B = - o(^/lO^ 

From the Report of Calibration for Platinum Resistance Thermometer 
Serial Wo. 1 781311 the following are the values for the constants. 

Rq , Resistance at the ice point 
= 25.5814 

Pundamental coefficient of coil 
0( = .003926833 

Delta in the Callendar formula 

1 .496176. 
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Actually the temperature t in °C at the outside of the tube 
protecting the platinum resister is given by 

t = t’ + M(t'). The small value M(t') is given by 


Ivi(t’) -.045 (i-^) 00* ' ^ ^419.58”'*^ ^ 


t' 


t' 


' 1 ). 


00 ' HI 9. 58 ^630.74' 

In our study this value ]l|(t') which serves to make the International 
practical scale of I 968 (IPTS-68) conform more closely to the 
thermodynamic scale" has been omitted. The omission was not 
intended but is due to overlooking. 


With the reported values of o( ^ Sand and neglecting 

M(t') a table has been generated on computer to give against 
temperature at an interval of .005°C, This table is later employed 
in directly reading off the temperature. It is intended to 
calibrate the thermocouple in the range 50°C to 180°C. A device 
that produces uniform temperature around the hot junction is the 
vapour jacket. This is a silvered, evacuated double walled unit 
(of dimensions I.D. = 40 mm. O.D. = 70 mm and length = 90 cm). 

The jacket is charged with 100 cc of pure liquid. The boiling 
flask is heated by means of a heating mantle. By such heating 
the liquid is brought in equilibrium with its own vapour. Under 
equilibrium conditions one can observe a constant drop rate from 
the tip of the thermov/ells for platinum resistance thermometer 
and the thermocouple. A steady state is reached when the 
recorder pen moves in a vertical line and there are very little 
or no fluctuations. Pressure is maintained constant. At this 
point the value of platinum resister is directly read on the 
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Recorder, emf in the thermocouple is read in mV from the potentio-i® 
meter. To alter the temperature conditions in the rapour jacket 
one has to change the pressure of the system by introducing 
Nitrogen via the toggle valves. 

Three liquids toluene, bromo benzene and naphthalene 
have bean used in series to cover the full temperature range 50°C 
to 180°C. Toluene has been used for range 50 to 110°C (corresponding 
to 100 mm to 760 mm pressure range). Bromobenzene has been used 
for 110°C to 155°C (i.e. 220 mm to 760 "mm) and Naphthalene for 
155®C to 180°C (corresponding to 1 50 * mm to 760 mm). Vacuum or 
low pressures have been maintained by plugging the vapour jacket at 
both ends vflth rubber stoppers which had been earlier soaked in 
kerosene to provide svjelling and 'therefore to get a close fit,. 
Apiezon Q was also used at both ends. Mercury was also used, at the 
bottom to provide a seal against air bubbles etc. The top rubber 
stopper has holes in which the thermowells for platinum resistance 
thermometer and thermocouple have been inserted. The holes were 
made by pressing under a drill machine immediately after immersing 
the stopper in liquid nitrogen. Since the platinum resistance 
thermometer was calibrated with a 7" depth of immiersion same depth 
has been kept under m.ercury in the thermowell. The thermocouple 
hot junction bead copper constantan has been formed by first 
baring, the wires from their sheaths, giving them three to four 
twists over each other with a nose pile- r, dipping in borax, and 
showing it to the flame till the bead is seen to be formed. The 
oxide that may have formed is gently rubbed off with hand . Similarly 
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the cold junction beads (Cu/Cu and const'antan/Cu ) are formed. 

• 

The enamel of the copper lead wire is removed by rubbing it off 
with sand paper. 

Trior to reading off the resistance values from the 

recorder, the pen of the chart is brought to zero by means of the 

zero rheostat. The zero point is also confirmed by reading when 

the thermometer is at the ice point. If the R value is not 

0 

exactly reproduced zero reheostat is adjusted to road the required 
value , 

Measurements of emf on the potentiometer are made after 
standardizing the current each time, '''ith reference to the' 
figure, below ? 


i 1- 



The operation of the potentiometer is as follows; ^/;lien the switch 
is placed on (1) then the current in the circuit is adjusted by 
means of rheostat C (Coarse and fine adjustments are possible) till 
there is no deflection in the galvanometer. Under this condition 
the current through the calibrated wire is "standardized". The 
switch is then Placed on (2) and the slide is adjusted until the 



101 


APPENDIX 2 

CALIBRAiTION OF BEPBACTOMEIER 

• , «r 

The anal 3 /sis of the vapour and liquid samples, collected 
after equilibrium has been reached, is done on an Abbe Refractometer 
v'ith the 4th decimal place being estimated by eye. The compo- 
sitions are then directly read off against the observed indices 
from the calibration chart. The calibration procedure is briefly 
described hero. 

Samples of knov/n composition are prepared. The weighings 
are done on an electronic balance which is accurate upto the fifth 
decimal place. The mole fractions are calculated by the following 
formulae 

W-| 

W2/m2 

^ w-|/m^ + W2/ni2 

i^Tiere w-j and W 2 are the weights of components 1 and 2 in the sample 
mixture and m-^ ,ra 2 are their respective molecular weights. About 
1 0 to 12 samples are prepared to cover the whole range of concen- 
tration. Prior to taking any reading with the Refractometcr one- 
must secure precise adjustment of the scale in its relation to the 
position of the total reflecting line. This is done by reproducing 
the index of the glass test piece. After proper cleaning of the 
prisms, a small drop of bromo naphthalene (used as contact liquid) 
is put on the lower prism and the test piece brought in contact 



102 


with it^with the engraved side up. The scale is set at the index 
val'^e engraved on the test piece. The iiluminant is swung in 
position directly in line with the prism surface, lamp shield ad- 
lusted to secure the hest contrast between the two halves of the 
field. The compensaLtor dial is turned until the dividing line is 
correctly compensated (red at one end and blue at the other). 

This position of tho compensator dial is then maintained for all 
readings . 

The prism is maintained at 30°C during all the readings by 
circulating water from a constant temperature bath. The bath 
temperature is controlled by on-off heaters and the required temperC' 
ture (30°C in our case) is set on a Beckman thermometer. The out- 
put of the bath is connected to the prism coupling in which the 
thermometer is mounted. The temperature was controlled to approxi- 
mately 1°C in order to obtain consistent and accurate results. 

Also it was necessary to keep the temperature constant for a 
reasonable period of time since the prism temperature always lags 
behind the bath temperature. 

Bor analysis of the samples, thorough cleaning, by an 
organic solvent, of the prism between samples is necessary. The 
liquid is introduced through a dropper and the prisms are then 
closed slowly squeezing out the excess. The scale is set at the 
approximate value expected and the borderline is brought upon 
the intersection of the cross hairs by means of coarse and fine 
hand controls. The index is then road by depressing the momentary 
contact switch, estimating the fourth place. All tho points for 
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calibration v'cre obtained in this manner. In addition the index 

of the pure components were also determined. 

The calibration data so obtained were then regressed by 

2 

applying least square technique on the polynomial a + bx + cx +... 
and terminating the degree of polynomial when the ratio of the 
variances which forms an P-distribut ion (experimental and control 
group) satis if ed the 95/^ confidence limit. 



APPENDIX 3 


EETAIES OE STMEABEIZATIOI 

In this chapter it is intended to follow up the develop- 
ment and evaluation of the apparatus, part by part, studying each 
minor detail both regarding construction and operation, so that 
we may have a fuller appreciation of our equilibrium study. 

The EquilibriumStill s 

The standardization of the- equilibrium still involved some 
ma;ior changes. With the first installation of the separation 
chamber and reboiler it was observed that the normal boiling point 
of pure carbon tetrachloride was not renroduced in the thermowell. 
The measured temperature, even after several hours of running and 
all precautions regarding insulation, was 73°C - about 5 to 6°C 
less than the true boiling point. Close examination revealed that 
there was ( i) no movement of vapour and liquid slugs in the Cottrell 
pump and (ii) no liquid seal around the thermowell in the separation 
chamber with the resultant escape of vapours into the liquid -side 
condenser. Seduction in diameter of the Cottrell pump from 10 ram 
I.E. to 7 mm I.D. brought about the much needed slug formation 
with resultant transfer of boiling liquid from the reboiler to 
the separation chamber. This brought about an increase in the 
temperature reading, for now the liquid around the thermowell 
was in rapid recirculation, and no longer formed a dead pocket. 
However, the vapour entrainment still remained a problem and this 
meant that the inclination of the thermowell had to be changed 
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so that liquid seal may be formed. Separation chamber was 
constructed a second time ?/ith the thermowell no?; making an angle 
of 50° to the horizontal line instead of the initial 10°. All 
the dimensions of the separation chamber were also scaled down by 
a factor of 8. The diameter of the orifice connecting the packed 
chamber to the thermowell was also reduced so that the hold-^up 
volume may reduce. 

le boiler size ?;as reduced to an effective volume of 160cc. 
The stop-cock for draining material from the reboiler was elimi- 
nated completely because the grease applied to it was being 
eaten a?7ay by the liquid, and without grease the stop-cork v/as 
leaking even at 760 mm Hg. A seal (4" high and parallel to the 
Cottrell pump) was provided which could be cut at the time of 
draining and the material taken out through a syringe. 

Under running condition of the equilibrium still the 
roboiler should be filled upto the neck (bottom tip of the 
Cottrell pump) so that the liquid may be drawn up the Cottrell 
pump in form of slugs. Before starting liquid level should be 
between 4" to 6" in the Cottrell pump. This required a total 
charge of 200 ’cc of the liquid mixture (40cc in excess of the 
actual volume of the re holler). The excess charge accounted 
for the hold-up in separation chamber and receiver -mixer under 
running condition. 

The performance of the Cottrell pump was measured by 
estimating the ratio L/V i.e. ratio of liquid to vapour recir- 
culation rate. And the Intensity of heating was measured by 
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e stin'^’t ing- the vapour ro circulation rate (counting the drop rate 
at tht- tip of the constriction provided in the vapour condenser)^ 
high vapour recirculation rate (proportionately controlled hy 
voltage) gave reproducible results for the pure components. But 
a beyond a certain value this caused flooding of the receiverr-mixer. 
This is the point of superheating. Flooding of the receiver- 
mixer should therefore he avoided at all costs. 

The equilibrium points for pure components were measured 
under the following conditions. 

L/Y V 

(ratio of liouid to (vapor recir- Heat supply Heat supply 
vanour recirculation cuiotion rate) to external to finger 

rate heater heater 

CyHi^ 3:1 200 drops/min. 20 volts 55 volts 

Cyll^^O 1 si 400 drops/rain. 25 volts 77 volts 

llhen data were collected starting with pure methylcyclohexane 
voltage was gradually increased and when it was collected from 
the pure oiethylcyclohexanol end voltage was gradually decreased. 

During the whole process of operation all the essential 
parts of the equilibrium still viz - separation chamber, receiver- 
mixer, re boiler, cottrell pump and the recycle lines are kept 
insulated by glass wool. Ideally the operation should b© performed 
in adiabatic surroundings. 

Single Stajfe Ebulliameter s • 

Like the equilibrium still the ebulliometer used for 
pressure measurement of the system has been lagged with glass wool. 
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Aboii*t 60 cc of pure benzene are necessary for slug fortnation 
in the Cottrell pump. The liquid is both filled in and drained 
out from the top of the condenser. This method eliminated the 
use of a stop cork for draining the liquid, it was necessary 
because the grease in the stop cork was being eaten away by 
benzene. In addition to the condenser an ice trap was provided 
at the top to prevent any escape of benzene vapours into the 
manifold. The ebulliometer is plugged to the manifold with a Bt2 
ground glass joint, 

A constant boiling point of benzene is reproduced for a 
condensation rate of 70 to 120 drops/min. The voltage was fixed 
at 1 9 volts to get a drop count of about 100 drops/min. The 
resistance of the heating element was 20 ohms. 

Differential Ebulliometer ; 

The whole differential ebulliometer has been lagged with 
.asbestos and magnesia powder mixture and on .top of it with glass 
Y?ool and asbestos tape. Amount of charge required for slug 
formation in the Cottrell pump is equal to 60 cc. A microheater 
has been provided near the top drop counter to boost the flow of 
vapors and prevent partial condensation occuring on the walls due 
to boat loss. 

The differential ebulliometer was standardized with respect 
to benzene. Both the thermowells indicated a temperature of 80,1 °C 
at 760 mm Hg when (i) main heater was at 30 volts and the micro- 
heater was at 6 volts and (ii) the drop counts were 30 drops/min. 
and 100 drops/min in the top and bottom drop counter respectively. 



